OS 62 


eo” May 1978 Volume 6 Number 5 


Popular Computing 


The maeceite for people who enjoy computing 


ae i. 
os 
6 
@ 9 
< 
qe 
gw 
sor 
wt 
41 we pt 
& 
aN) 10; 
5 
3 
@ 


PC62--2 


A racing driver travels a circular course, taking 
the following steps: 


1/2 a revolution, then 

2/3 revolution, 

3/4 revolution, 

4/5 revolution, 

5/6 revolution,... and so on. 


At the end of each step he waves to his girl, at START. 
As shown by the figure, at the 10th step he is close to 


her (within 7.156 degrees). At the 29th step, he is 
even closer. Eventually, he can get as close as he 
pleases. 


The Problem is, at what step will he be within some 
arbitrary distance of START, say one degree? 


Contributing Editor Thomas R. Parkin writes: 


"The problem is to sum the series 1/2 + 2/3 + 3/4 +... 
for some number of terms until the sum is within 1/360 
of an integral number of revolutions. This can be done 
by evaluating a cumulative sum of n/ntil, starting at 
n= 1, and continuing until the fraction is close enough 
to an integer. 


We could also note that this sum can be expressed as: 
eet ett 1/3) + (1.- 174) + ... = N= Din) 
The second term of this sum is the sum of the harmonic 


series and unfortunately there is no simple formula for 
the sum of N terms of the harmonic series. 


PROBLEM 2 2 9 


The following program, written in the JOSS language, 
will handle the problem as given. Everything in the 
program is self-explanatory except the function in step 
it tis eaGoyy This stands for ‘integer part." 


Continued on page 17... 
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Contest 15 Result 


Contest 15 (Reverse Knockout ) in our December issue 
brought entries from all over the world. 


The Problem was this (see Figure 1): Start with all 
the natural numbers from 3 up. Reverse all sets of 3. 
This brings 5 to the head of the stream, which dictates 
knocking out the 5th number after the 5, which is 6; the 
5 replaces the number knocked out. The process now 
repeats. The leading number is 4, so all sets of 4 are 
reversed, bringing 7 to the front, and the 7 knocks out the 
Die The list of numbers that are knocked out (those 
circled in the diagram) is to be extended. 


When the problem was constructed, it was our belief 
that the only possible solution would involve actual 
manipulation of huge blocks of numbers; that is, a brute- 


force approach. This would tend to restrict entries for 
the contest to those who have access to large and fast 
machines. 


Associate Editor David Babcock wrote a program in 
assembly language that produced a list of the first 48 
numbers, using this crowbar technique; the run took 36 
seconds of CPU time. A longer list, using the same 
& approach, would require much storage swapping and overlays 
and enormous amounts of machine time. 


We were pleasantly astonished, then, to find that 
most entries involved a combined analytic and computer 
attack on the problem. The material submitted by our 
contest winner: 


John D. Beeby 
Milibrae, California 


is an outstanding piece of computing, and we are pleased 
to display most of it verbatim. 


But even the entries that chose to manipulate lists 
of numbers surprised us with their ingenuity, elegance, 


and sophistication. The best example is from Joan Farmer, 
of the Polytechnic of North London (the winner of our 
eleventh contest). It is worthwhile to study her program 


and her analysis: 


"The best solution to this problem is an analytical one, 

Having failed, as I did, to find such a solution, one is faced 
& with the impossibility--common to this type of problem--of 
packing an unbounded list of numbers into a finite machine. 
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"My first approach was based on the fact that we are 
continually being told that hardware is fast becoming so 
cheap that we can more or less disregard it, while programming 
costs are continuing to soar. So I decided to forget about 
economies of space and time and concentrate on producing 
a simple, clearly constructed program that would mirror the 
inherent features of the problem and be easy to understand 
and maintain (even if it did need a computer with an out- 
rageously large memory to produce a reasonably long list of 
results). 


"The problem revolves ‘round two lists--a list of the 
natural numbers and a list of "leaders"--both of which are 
continuously growing and need to be permanently remembered. 
The list of numbers contains the integers in their current 
state of reversal, and each leader remembers the span of its 
reversais and how far through the numbers these reversals 
have been made. No leader may, of course, make a reversal 
beyond the point reached by its predecessor. Thus, a 
natural need for recursion arises. 


"I decided to hold leaders and numbers as linked lists, 
the links of the number list being altered to record 
reversals rather than swapping actual values, and both lists 


being held in internal memory. The program proved to be 
extremely simple. It was written in Algol 68 which provides 
the necessary list processing and recursion. Each new 


leader reverses his first span of numbers and then continues 
reversing until his candidate for knocking out is wifhin 
range. If necessary, he calls on the previous leader to do 
sufficient reversals to cover the required range. If the 
recursive calls for reversal get back to the original 
leader, he extends the list of numbers by adding on reversed 
groups of 3. 


"Only four procedures were required: 


extend--to add new numbers to the end of the 
number list; 

reverse--to provide the recursive reversal 
mechanism; 

kKnockout--to pick off the next number for printing 
and replace it by the current first 
number; 

newleader--to introduce the next leader to the 
head of the list. 


"There was only one snag with this program. it 
didn't have the postulated cheap but extensive hardware to 
run it on and so it used up all the available memory after 
printing only 41 numbers. I decided to submit my solution 
and think again!" 


Joan Farmer's 
Algol 68 
program 


knockout pe Dec 77 


= struct (int n, ref number next); 
= os 


= struct (int val, c¢ length of exchange 
start, ¢ ordinal of next chain 
ref number last, c¢ end of last chain 
ref leader next); ¢ previous leader 
eNO Teerelese joni ) yonlale 
ref number niln = nil; 
ref number firstn t= number := (2, niln); 
ref number lastn := firstn; 
ref leader first] := leader := (3, 3, firstn, nil); 
nt a i= 3, Zz := 2; c ordinals of start and end of integers 


proc extend = void: 
e add 3 more numbers to the number list in reverse order c 


begin 
for j from z plus 3 by -1 to z- 2 
do begin 
next of lastn := number := (j, niln); 
lastn := next of lastn 
end 


end; c¢ end of extend c 


proc reverse = (ref ref leader 1) void: 
e reverse the order of the next set of numbers in 1's chain 


begin 
ref number 1f, n, p, q; P3 


iieenextuor lis nil 

Then (extend; Last of 1 := lastn) 

else while (start of 1+ val of 1) >= start of (next of 1) 
do reverse (next of 1); oF 
If := last of 1; 
pi> next of lf; q := next of'p; mr := next of q; 


—_ 


7= p; c¢ hang on to the new 'last!' c 


ago begin 
next of q := p; c¢ reverse a pointer c 
Pp 3= G5 q :=r; ¢ and Cc 
roi= next of r © move on one. c 
end ; 
next of q := p; ¢ tidy up ends @ 
next of n := 7; 


last of 1 :=n 
fi; 
Start of 1 plus val of 1 
end; © end of reverse ¢c 


jolelelo 


Jo 
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proc knockout = (ref leader 1) void: 
@ print next candidate and change its value @ 


begin 


ref number rn := firstn; 
to n of firstn do rn := next of rn; 
print (n of rn); 
if charsleft(standout) 10 then print (newline) fi; 
nof rn :=n of firstn 
end; ¢ end of knockout c 


proc newleader = vold: 

c¢ add a leader to the head of the list of leaders c 

begin 
@irecie:= leader := (n of next of firstn, a plus 1, firstn, firstl); 
firstn := next of firstn c drop the first number c 

end; ¢ end of newleader c 


¢ the main program starts here c 


print (("Popular Computing contest 15(PC57-2, Dec 77)", 
newline, newline) ); 
reverse (firstl); 
firstn := next of firstn; 
do begin 
reverse (firstl); 
while a +n of firstn >= start of first] do reverse (firstl); 
knockout (firstl); 
newleader 
end 
end 
finish 


But the brute-force approach, however cleverly 
implemented, can never extend the knockout list very 
far; after relatively few stages, the procedure would call 
for reversing blocks of millions of numbers and even with 
extensive (virtual) storage, the CPU time would become 
excessive. So the analytic approach must take over. 
What follows, then, is Mr. Beeby's solution. 


In my solution the diagram of Problem 212 (PC57-1) 
is expanded, figuratively speaking, one "cycle" (that is, 
one pair of lines) at a time down the page and as far to 
the right as necessary. Actually, all of the numbers in 
this vast hypothetical structure are ignored except for 
three key elements in each cycle. One of these is the 
extraction or output item. The others (described in 
Figure 2 as “starter” and "pointer") are stored for use 
by the program, Hence the modest storage requirement is 


of three key elements is determined is explained in detail 
in the discussion which follows. 


two items per output line. The method by which each set » 
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I validated my analysis with a BASIC program, which 
I converted to Fortran in order to use double precision, 
produce a longer output list, and check the execution 
speed on a big machine. The Fortran H-Extended venture 
gave a still longer output list and that was its sole 
contribution to the project. I am indebted to the 
College of San Mateo and Stanford University for the 
privilege of using their computer facilities. 


10 REM * PROBLEM 212 JOHN D. BEEBY 


20 DIM S(100),P(100) 
30 N=O 

4o eRe 

50 P(1)=5 

60 N=Nt+1 

70 V=N+P(N) 


So) eapacge ON oe 
100 IF N=100 THEN 330 
110 V=Ntl 

120 GSUB 190 

130 S(N+1)=x 


John Beeby's 140. V=N+X 
BASIC program 150 G@SUB 190 
160 P(Nt1)=x 
170 GTS 60 
180 REM * "TRACE" SUBRZUTINE 
190 I=N 


200 C=V-I+tl 
210 Q=C-INT(C/S(I))*S(T) 
220 IF Q >= 1 THEN 240 

230 Q=S(I) 

240 V=V+S(I)-@-Q+t1 

250 IF I <> 1 THEN 280 

260 X=V+2 

270 RETURN 

280 IF V <> I-1+P(I-1) THEN 310 
290 X=P(I-1) 

300 RETURN 

310 I=I-1 

320 GPTP 200 

330 END 


"Trace" Rationale 


If we define the location of some element as Column V, 
Cycle N and it is in the Qth position of a reversal field 
of width S(N), then it moved 


S(N) - (2Q - 1) 
positions as a result of the Cycle N reversal, I call 
the above expression the "swing." When evaluated, it 
indicates how far to the left (-) or right (+) we have to 
move in order to arrive at the previous location of the 
element. 


D 
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If an element is in the extreme right hand position of 
its reversal field, Q is equal to S(N) and the "swing" is 
equal to -S(N)+1. That means that the previous location 
was S(N)-1 positions to the left, and the element was in 
the extreme left hand position before the reversal. 


If an element is in the extreme left hand position @ 
of its reversal field, Q is equal to 1 and the "swing" is 
equal to S(N)-1. We find the previous location by moving 
S(N)-1 positions to the right which places the element in 
the extreme right hand position before the reversal. 


Of course, if S(N) is odd and Q happens to be equal 
to (S(N)+1)/2, the “swing” will be zero, placing the element 
in the same position before and after the reversal. 


In general, 


V(N-1) = V(N) + S(N) - (20-1). * 


By this reasoning we can trace upward through the 
chart to reach the top line, but there is one special case 
to consider at each cycle on the way. The present value 
in any former extraction site did not arrive there in the 
usual manner, through a series of reversals. Instead, 
this value moved in from the “pointer” logation of the 
previous cycle. Thus, when the trace routine lands on an 
extraction site, the appropriate value can be recalled 
from array P and the trace is done. 


*In the programs, this translates to V=V+S(I)-Q-Qtl. @ 


The following names appear in the programs and in 
Figures 2 and 3. 


N is the cycle number. 

I is a stand-in for N. It is set equal to N at the 
beginning of the Trace routine and then decremented 
during the operation of the routine. 

S is an array dimensioned for storing the value of S(N) 
corresponding to each output item. S(N) is the 
"starter," the first element in the first line of 
cycle N. (In the original problem specification this 


was called the "leader" and designated K.) The value 
of S(N) defines the width of the reversal field 
in cycle N. 


P is an array dimensioned for storing the value of P(N) 

corresponding to each output item. P(N) is the 

"pointer," the first element in the second line of 

cycle N. It "points" to the extraction site. 

is the column number, referring to the entire chart. [ @ 


ag 


is the column number within a specific cycle. Cycle N 
begins in column V, where V is equal to N and C is equal 
to 1. C is always equal to V-Ntl (or V-I+1) so it is 
just an tntermediate to simplify another expression. 


Q is the position of an element within its reversal field. 
Q refers to the second line of the cycle, after the 
reversal has taken place. Q is equal to C mod S(N), 
with the proviso that if Q = 0, then the value of se) 
is assigned to Q. 

X is the value returned by the Trace routine in the BASIC 
program and the flowchart. In the Fortran program as 
well as the BASIC program it designates the extraction 
value in the WRITE statement. 


The following additional names appear in the Fortran 
program. They were included to minimize conversions 
and avoid mixed mode expressions. 


is the Real*8 equivalent of N, 
is the Real*8 equivalent of I. 
is a constant, 1.0D0 . 


GAaK 


Comments on the Programs 


The BASIC program was written in Hewlett-Packard 
Time-Shared BASIC/2000, Level E, and run on the HP-2000 
at the College of San Mateo. The output list was 
limited to 100 items to keep within the number size of 
the BASIC system. Execution time was approximately 
3.5 minutes. The main program or the trace loop or both 
could be written as FOR...NEXT loops. Here and also in 
the Fortran program I preferred to use explicitly coded 
loops. This allowed the BASIC program to conform exactly 
to the flowchart and also avoided some minor loose ends 
that seemed to be inherent in the FOR...NEXT approach. 


The Fortran program was written initially for the 
IBM Fortran G compiler and run as a compile~and-go batch 
job on the IBM System 370/168 at Stanford University. The 
program had a net count of 62 source statements. With 
double precision, accuracy problems began to appear around 
the 250th line of output. The program executed in about 
1.95 seconds of CPU time. 


The final program was run with the Fortran H-Extended 
compiler, with the necessary program changes to specify 
Real*16 for the variables; the compiler optimizing options 
were not used. The run for 350 lines of output (notice 
that the solution got to 23-digit numbers) took 36.39 
seconds of CPU time. 


The problem posed in Contest 15 is intrinsically 
useless. It seems to be in the nature of the world that 
it is only such problems, when they are sufficiently 
challenging, that attract the best minds and the most 
effort. 

The 350 results that Mr. Beeby obtained, then, are 
not in themselves of any value, but we wish to present them 
for the record. 
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101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
abil} 
114 
aS) 
116 
LAG, 
118 
119 
120 
U2 
122 
123 
124 
125 
126 
ley 
128 
129 
130 
IS} 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 


386162 
1329220 
797702 
1329212 
821084 
895228 
3188714 
1851050 
14832 
31322 
458010 
3688224 
116162 
3166184 
2227184 
1121388 
5578302 
8772416 
8247148 
8772622 
8247218 
8772686 
8086144 
3251348 
4060962 
3251062 
3556682 
14894278 


20523970 


21792708 
43116530 
45675326 
71801724 
32044684 
51378470 
89684842 
1023048 30 
342784454 
319528456 

6615548 
11450782 

6614898 
11406830 
16197704 
74724060 
150807696 
217464252 
135382098 

7301960 
43535338 


15: 
Loe 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
alyjal 
1/2 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 


62096664 
506484280 
204962322 
163462412 
398156956 
871279960 
540025168 
575884716 

1047167076 
2178029746 
4231263224 
4117932012 
14586 30740 
1195303030 
2200164448 
105764 3640 
851838540 
3119982314 
4055121086 
5087750994 
6560447254 
17055331578 
18983964 368 
38249123480 
64571911298 
28096848050 
34244417480 
4043762140 
26504033170 
37033757256 
86182666702 
24865251648 
59674595214 
46679361438 
156727865474 
195380136622 
183843194272 
234513417442 
231375004856 
2447609219 36 
83566681700 
288741029076 
658615695470 
97 3492533902 


1197324615328 
1160943169758 


939909285014 
816942516594 


1700983282472 


873368821832 


The 350 results to the Reverse/Knockout 
problem obtained by John D. Beeby. — 


» 


234 
235 


1647940441800 
629077431134 
1515431095038 
1094423372758 
439697 24 3988 
448533959454 
2765040074696 
1588989792664 
2930116079276 
2724480645010 
2170128551744 
795478898464 
1618908577662 
6899877 34 37 38 
9227980445560 
6899877 3460/44 
9227979400318 
155495967454 34 
6169390800316 
13054046894786 
39167901725720 
13054046916796 
39167901696180 


746329181357 34. 


37597628692670 
10033206912682 
56444 394006990 
10033206912568 
56445261921258 
14409638794 316 
413905390285006 
5486961 3815198 
149037264019094 
4454 3859198258 
15618 3522614262 
3346764 31756840 
4970023897 580 36 
80 3243570041582 
1367 520555923724 
1689166854660282 
2764.67 6536466134 
2266659858104154 
915251093748280 
12994484256687 30 
3986 391157615022 
1299448425668858 
3986 39115487 3888 
8717089 329013488 
867219 3321601416 
8756264978594208 


94. 32599231110202 
220088 29507645552 
13942311078213494 

5713338879885968 

502476220746 3218 
15328397915860412 
21145341980229790 
344 36 364742009660 

870516 323968750 

3393279985709486 
6560274044 3142098 

3393279985711506 
6560274044 3104408 
33131777600241056 

29190078 3555120 

7155 30633647864 
42037420106576402 
26448205 349489448 
6897985075747 3234 
465621 36328422 358 
1084 34261541415110 
465621 36328422 322 
1084 342615414052 34 

1722328149420166 
31928502686558550 
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479218831512570 3600100 
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Set m= 1. 

Set t = 0. 

Set n = ntl 

Set m = ml. 
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UARES WITH RENDOM SCALRACTERS 


The Affectionate Racer 
Continued from page 2. 


Further calculations revealed some of 
the critical values for the parameter, c, as 


shown in this table. 


The next critical 


value would require very high precision 


arithmetic. 


t 


1.1666666 
1.9166666 
7.0710318 
79801227 
26 .0050129 
TT 997932 
220 .000039 
4540 .00001 


| 


. 1666666 
-0833333 
.0710318 
-0198773 
-0050129 
.002068 
-000039 
-0000095 
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El 


We begin here a new feature: homework assignments for 
computing students. Each assignment has been debugged 

and time-tested through many semesters of use. Most 

of the loopholes have been plugged, so that a student r 
should know what to do from the directions given here. 
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Differencing 


Given a deck of cards. Each card bears one number-- 
a single value of some function. 


Draw_a flowchart for the following logic: after 
reading each card, print on one line the functional value 
and the first, second, and third differences, 


The length of the deck is unknown, but you can assume S 
(since we want to go to third differences) that it is at 
least 10 cards long. The end of the deck is detected by 
the simplest of end-of-file procedures; namely, that there 
are no more cards to be read. 


For an elegant solution to this problem, you could 
spend a lot of time worrying about not developing or printing 
various differences on the first three cards, since they will 
be garbage anyway. You are to ignore this fact; the garbage 
can be crossed off the listing with a pen. If that sounds 
like sloppy computing, keep in mind the dictum "First make it 
work--then make it pretty." We are concerned here with the 
first part of the dictum. For now, we want only to get the 
differencing logic expressed. 


When you have completed your flowchart (which should be 
fairly simple), consider the following points: 


Have you made any assumptions about the data, such as 
its order, or whether the numbers are all positive? (There 
should be no assumptions.) Does your flowchart produce 
differences, as asked for, or have you produced the absolute 
values of the differences? How much storage would your 
logic require? Will it still work if the data deck is 


100,000 cards long? Have you identified all your variables? 
Does your flowchart require a personal lecture to go with it, 
or could anyone understand it? This is intended to be a 
simple, straightforward task--don't make it into something 
grandiose. 


ODD & EVEN 


An integer is stored at address X. we have come 
to the point in a problem solution where our course of 
action depends on whether the number at X is odd or even; 
that is, to the situation expressed on a flowchart as: 


The direct and natural thing to do is divide X by 
2 and test the remainder of that division for zero (and 
if zero, X is even, and we go to reference 3). The 
trouble with that is 


1. division is a costly operation in terms of 
CPU time, and 


2. for most coding languages, it is messy 
to extract the remainder of a division. 


There is, of course, the MOD function in many 
languages which does what we want, but how does it get 
implemented at the machine language level? 


Similarly, some languages have op-codes for or) 

JUMP ON ODD but again--how do they work? “ 
a 

Try another tack. Perform an AND between word X P| 

and the number one. This operation will mask in the unit's re} 
& 

Qa 


bit of X in the accumulator, so that a JUMP ON ZERO will be 


taken if X is even. 
There must be many other ways... [dal 
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GENERATIONS 
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Vacuum tube machines. Early models one-of-a-kind. 
First mass produced machine the Univac I. First, 
drum and CRT storage; cores on later models. 
Included IBM 650 and 704. Decimal machines for 
business; binary for scientific. Reliability 


1950-57 relatively low. 


Transistorized machines. Included IBM 7090, 1620, 
and 1401, IGP-30 and G-15 dominated the small 
machine market. Still had decimal equipment for 
business purposes. 


58-64 
Hl e 


Byte logic, and 1/2" chip circuitry, typified by the 
System 360. Now all machines are binary. The 
16-bit minis began around 1968, 


64-70 


Floppy disk control, typified by the System 370; for 
the first time, identical machines were operating in 
the field. 8-bit micros began in 1973, and the 
personal and hobby machines started in 1975. Rise 
of virtual storage. 


The watchword is firmware, as more and more of the 
burden of overhead is shifted from software to a 
hardware. 50 nanosecond complete cycle times, 

and going down. Still wedded to magnetic cores, 

at least at the start of this generation. Power- 

ful systems being built around microprocessors. 


